Dynamic structure function of a cold Fermi gas at unitarity 
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We present a theoretical study of the dynamic structure function of a resonantly interacting two- 
component Fermi gas at zero temperature. Our approach is based on dynamic many-body theory 
able to describe excitations in strongly correlated Fermi systems. The fixed-node diffusion Monte 
Carlo method is used to produce the ground-state correlation functions which are used as an input 
for the excitation theory. Our approach reproduces recent Bragg scattering data in both the density 
and the spin channel. In the BCS regime, the response is close to that of the ideal Fermi gas. 
On the EEC side, the Bose peak associated with the formation of dimers dominates the density 
channel of the dynamic response. When the fraction of dimers is large our theory departs from the 
experimental data, mainly in the spin channel. 
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An impressive advance in realizing and controlling ul- 
tracold Fermi gases has permitted to study physical phe- 
nomena whose appearance was previously only a mat- 
ter of speculation |l|. One of the major advances has 
been the first clean observation of a crossover between 
two limits of Fermi matter, the Bose-Einstein conden- 
sate (BEC) regime composed by Fermi molecules and the 
Bardeen-Cooper-Schrieffcr (BCS) superfluid gas where 
Cooper pairs are formed [^-Isl- 

In the dilute limit, interactions between atoms can be 
described by a single parameter, namely the s-wave scat- 
tering length a. The system evolves from a BCS regime 
with no bound state (a < 0) to a molecular one (BEC) 
with a two-body bound state (a > 0), crossing a singular 
point where |a| — > cx) corresponding to a Fano-Feshbach 
resonance [6|-|8[. This special point is referred to as the 
unitary limit and is expected to show a number of uni- 
versal properties. Indeed, in this regime the only relevant 
length and energy scales are the inverse Fermi momen- 
tum kp and the Fermi energy Ep, respectively. 

Nowadays, the unitary Fermi gas is routinely realized 
in experiments. An acccurate theoretical description of 
the unitary limit requires the use of advanced methods 
of many-body theory that can deal with strongly inter- 
acting systems. In this situation, quantum Monte Carlo 
methods have proven to be very useful for the calculation 
of the finite- and zero-temperature equation of state (a-Q 
which is found to be in good overall agreement with ex- 
perimental data lOljlUl- Also, it was possible to predict 
the effect of the unbalanced population [12| and different 
masses [l3'|. More recently, static responses such as the 
spin susceptibility and spin diffusivity |14l . |l5| have been 
measured and also studied theoretically using quantum 
Monte Carlo (QMC) methods [16| and a diagrammatic 
approach [17| . 

Recently, Hoinka et al. 18|, ll9| have measured the dy- 
namic spin and density responses of a Fermi gas at uni- 



tarity and in the BEC regime. In order to measure both 
responses, they used two-photon Bragg scattering, acting 
on ^Li at very low temperature (~ 0.06Ef at unitarity), 
and a proper choice of the laser detuning. The density 
response at unitarity is significantly different from that 
of an ideal Fermi gas, showing a clear peak which corre- 
sponds to the formation of dimers composed of spin-up 
and spin-down particles. This signature of dimer for- 
mation is not visible in the spin channel, nevertheless 
our results are quite different from the ideal Fermi gas 
response. This provides evidence for the importance of 
correlations in the unitary gas despite its diluteness. 

Some of the features shown by the dynamic response 
function at unitarity were previously determined in the- 
oretical work: at T = using a dynamic self-consistent 
mean-field approach based on BCS theory [20| , at T > 
using virial expansions [21]. However, the dynamic prop- 
erties along the BCS-BEC crossover are not completely 
understood due to the need of a full many-body the- 
ory and the impossibility to use QMC methods directly. 
The goal of our study is to pursue further a theory for 
the dynamic structure function, combining the virtues of 
Monte Carlo methods and modern diagrammatic many- 
body theory. 

In the present work, we study the density and spin re- 
sponses of the two-component Fermi gas at unitarity and 
in characteristic points on both the BEC and BCS sides 
of the crossover. We use a fully microscopic approach uti- 
lizing correlated basis function (CBF) theory. The input 
ground-state structure functions are obtained from diffu- 
sion Monte Carlo simulations within the fixed-node (FN- 
DMC) approximation [6|. In its most advanced form [22], 
the dynamic many-body theory (DMBT) used here is the 
fermion version of the CBF-Brillouin-Wigner perturba- 
tion theory previously developed for bosons [231 ]. The 
power of DMBT for fermions has been recently demon- 
strated by the prediction of a stable roton excitation in 



two-dimensional ^He, which has independently been con- 
firmed by inelastic neutron scattering experiments [2J|. 

Inputs for the calculation of the dynamic response are 
the static structure factors, <S'-|-|-(fc) and Sf^{k), and the 
fraction of dimers along the BCS-BEC crossover. We 
obtain a complete description of the density and spin 
dynamic responses in the momentum-energy plane. In 
the unitary limit, both theoretical responses are in sat- 
isfactory agreement with the recent experimental Bragg 
scattering measurements [18| . 

In weakly interacting systems, time-dependent 
Hartree-Fock (TDHF) theory [2|] is a well-established 
method for capturing the dynamics. One assumes that 
the system is subjected to the perturbing Hamiltonian 



SH,^t{t)= / dV/"^(r)/iext(r,i) 



(1) 



the superscripts p and a stand for density and spin ex- 



citations, respectively. p^''^'^\r) is the (spin-) density op- 
erator, and /iext(r, i) is a weak, local external field. The 
wave function resulting from such a perturbation is as- 
sumed to be of the form 



\^) = e^/^^pl.'^p'^Wala^^^^^ ^ 



(2) 



where \(j)o) is the ground state of a non-interacting sys- 
tem with the same density. The particle ("p") and hole 
( "/i" ) labels run over spatial quantum numbers and spin 
degrees of freedom. 

The amplitudes Cph{t) are Fourier decomposed 

<^ph{t) = 4t^(^)e-^("+"'/'"')* + c'^^l^*{u:)e<^+'^"'^' . (3) 

Equations of motion are then derived from a least action 
principle 26|, [27[ . They are given by 



p'h' p'h' 

E y^;vP'>^h-c% {lu) =2 / d^pl^l;:^ (r)/iext(r, c^), (5) 



p' h' 



where /icxt(r,a;) is the Fourier component of the external 
field. 

From the solutions cL (cj) we can calculate the (spin-) 
density response function 



and the (spin-) density dynamic structure function 
S'(''/'^)(r,r',w) = -ilmx^''/'''(r,r',w) 



I 

the Hartree-Fock single-particle energies. Using the CBF 
method, the descri ptio n can be extended to strongly in- 
teracting systems [28i, '29'. The basic outcome of the 
rather intricate diagrammatic analysis is that one arrives 
at the same equations of motion. CBF theory provides 
a method for computing effective interactions from the 
bare interaction. These interactions are normally differ- 
ent in the direct and the exchange channels. 

If one neglects exchange effects, one can construct Vd 
from the static structure function by demanding that Vd 
(') reproduces S^''^q) through the uj'^ sum rule 



(6) 



in coordinate space. In a translationally invariant geom- 
etry, their momentum space representations, x^''''^' {q, t^) 
and S^P''^' {q, w), are functions of the momentum transfer 

Key ingredients of the equations of motion are the ma- 
trix elements of the interparticle interaction 

Vph-p'h' = {ph'\ Vd IV) - {ph'\ T/cx \p'h) (8) 

Vpp'hh';0 = {pp\ Vd \hh') - {pp\ Vex \h'h) (9) 

and the single-particle spectrum tph = ep — e^- Gener- 
ally, the potentials Vd and T4x are effective interactions 
in the direct and exchange channel, respectively. In con- 
ventional TDHF theory, they are the same and just the 
bare interaction between the particles. The Cp, eh are 



S'^p\q)^ / d{huj)S'^P\q,uj) 



(10) 



A useful analytic relationship which is accurate within 
one to two percent can be derived using a collective ap- 
proximation for the Lindhard function. 
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(11) 



where Sp {q) is the static structure factor of the free Fermi 
gas. 

For spin-independent interactions, the spin-density 
fluctuations depend on Vcx{q) but not on Vd{q). Thus one 
might be tempted to construct Vcx{q) from S^'^'{q) in a 



similar fashion as Vd{q) is derived from S^^^^q). Unfortu- 
nately, this procedure has numerical difficulties because 
the resulting problem is very poorly conditioned. We 
therefore utilize an expression obtained from the Fermi- 
Hypernetted-Chain (FHNC) analysis which suggests, in 
its simplest form [30|, a local exchange interaction 



14x(g) 



?iV S(P\q)-SF{q) 



2m 



SFiqY 



(12) 



The same effective interaction defines the single-particle 
spectrum 



2 1,2 



Efc 



h'k 



2m vN 



;:!^E^e.(ik-h| 



(13) 



where v is the spin degeneracy and N the number of par- 
ticles. A consistent treatment of the exchange potential 
and the single particle spectrum is needed to satisfy the 
ijj'^ sum rule. 

The calculation of the potentials Vd{q) ([TT|) and 
Kx(9) ([T2|) requires knowledge of no more than the den- 
sity structure function of the interacting system, S^''\q). 
For that, we use ground-state results obtained using the 
FN-DMC method [6|, which has proven its accuracy in 
the description of the BCS-BEC crossover. 

We have calculated the dynamic structure function in 
the density and spin channels for l/{kFa) = — 1 (BCS 
like), l/ikpa) = (unitarity) and l/ikpa) = 1 (EEC 
like) including exchange terms and non-local, energy- 
dependent CBF corrections that are not spelled out ex- 
plicitly, see Ref. [29|. Paying attention to the fact that, 
with decreasing l/(k-pa), the system favors the creation 
of bound dimers, we consider a non-interacting mixture 
of dimers ("bosons" with 2m) and atoms (fermions). The 
relative concentrations by which we weight the two con- 
tributions to the dynamic structure function in the den- 
sity channel are obtained from our FN-DMC ground- 
state calculations. In particular, the concentration of 
dimers as a function of the scattering length is estimated 
assuming that all pairs are in the zero-momentum state. 
According to this criterion, the concentration of dimers is 
100% deep inside the BEG regime and decreases to zero 
quickly after crossing the unitary limit where it amounts 
^ 50% [31| . Our model cannot predict the width of the 



"Bose" peak, hence it has been fitted to reproduce the 
width of the experimental data subject to the constraint 
that its contribution to the a;° sum rule is constant. The 
spin channel dynamics on the other hand is not affected 
by the presence of dimers because they do not contribute 
to spin-density fiuctuations. 

In Fig. [H we show results of S^P^I^k^ui) and S'^''\k,uj) 
for the wave vector q = 4.5 kp- We have selected this 
particular q value because it is the only momentum mea- 
sured by Bragg scattering [l8i|. In agreement with the 
experimental data, the density channel of the dynamic 
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FIG. 1. (Color online) The dynamic structure function in 
the density (top) and spin (bottom) channels at unitarity 
(l/(kFa) — 0) and at a momentum q = i.bkp (solid hue). 
The points are the experimental data of Ref. [13|. Dashed 
line stands for the non-interacting response. 



structure function shows a clear "Bose" peak centered at 
energy h'^k'^/Am coming from the scattering off the t-4- 
pairs. At higher energies, contribution of single Fermi 
atoms appear as a broader response. The importance of 
correlations at unitarity is best demonstrated by com- 
paring the response coming from the full theory with the 
non-interacting Lindhard function, also shown in Fig. [1] 
Our result for the spin channel at unitarity and the same 
q value as in the experiment is also shown in Fig. 1 
(bottom panel). As mentioned above, dimers do not con- 
tribute to the spin channel and only atomic contributions 
appear. As in the density channel, the use of many-body 
theory is crucial to describe the experimental data, note 
that the response of a non-interacting system has a com- 
pletely different structure and differs substantially from 
what was observed in the experiment. 

The full momentum-energy dependence of the density 
and spin response functions is shown as density plots in 
Fig. [2] At low momenta, k < Akp , the dimer peak lies in- 
side the particle-hole (p — h) band and has a large overlap 
with the single-atom scattering contribution. At higher 
momenta, the "Bose" peak is progressively decoupled of 
Fermi excitations with a position given by half the recoil 
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FIG. 2. (Color online) The dynamic structure function in 
the density (top) and spin (bottom) channels at unitarity 
(l/(kFa) = 0) as a function of the momentum transfer k and 
energy transfer fuu in units of the Fermi energy Ep . The solid 
lines show the limits of the particle-hole band. The vertical 
dashed line corresponds to the momentum transfer at which 
the experiment of Ref. Q^ was performed. 



FIG. 3. (Color online) The dynamic structure function in 
the density (top) and spin (bottom) channels at unitarity 
(l/(kFa) = 0), BEC (l/(kFa) = 1), and BCS (l/(kFa) = 
— 1). The lines are results from dynamic many-body theory 
and points are experimental data of Ref. [l^. All data corre- 
sponds to g = 4.5fcF. 



energy. The strength of the spin response is localized in- 
side the p—h band but interactions make the peak slightly 
asymmetric and shifted to higher energy compared with 
the non-interacting response; this effect is washed out at 
momenta k > ikp. 

Bragg scattering 18|, ll9| was also applied to a point in 
the phase diagram located in the BEC side, in particu- 
lar at 1/ (kpo) — 1, whereas no experimental data exists 
in the BCS side. We have applied dynamic many-body 
theory to the experimental BEC point and to the sym- 
metric one in the BCS side, with l/(fci?a) — —1, and at 
a momentum transfer q = A.bkp. The response in the 
BCS regime is rather close to that of the free Fermi gas, 
with only an incipient Bose peak from the dimer contri- 
bution due to its very low concentration (< 10 %). On 
the contrary, the dimer contribution dominates the den- 
sity channel in the BEC side since the concentration of 
dimers is ~ 90 %. Our theory, which produces results 
in good agreement with experiment in the unitary limit, 
seems to worsen in the BEC regime. This is seen most 
clearly in the spin response at l/{kpa) = 1 where our 



result does not show the broad response obtained from 
Bragg scattering. 

To conclude, motivated by the first measurements of 
the dynamic structure function of a two-component res- 
onantly interacting Fermi mixture, we have developed a 
theory able to predict its dynamic response. To this end, 
dynamic many-body theory is used together with a FN- 
DMC input for the static structure factor. From a direct 
comparison with experimental data [18j . we show that 
the ideal Fermi gas model is not sufficient to reproduce 
correctly the dynamic response functions. We have for- 
mulated a description in terms of a mixture of dimers 
(bosons) and atoms (fermions). This simple model per- 
mits to capture correctly the complicated structure of the 
density response at the unitary limit. Our theory departs 
from the experimental data, mainly in the spin channel, 
when the fraction of dimers is large. Finally, we provide 
information on the dynamic response for different (pre- 
viously not measured) values of the momentum in the 
unitary limit and on the BCS (not measured) side. 
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